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SUMMARY 


A parameter estimation technique has been used to analyze the August 1975 
Nimbus 6 Earth radiation budget (ERB) data to demonstrate the concept of decon- 
volution. The longwave radiation field at the top of the atmosphere is defined 
from satellite data by a fifth degree and fifth order spherical harmonic repre- 
sentation. The variations of the major features of the radiation field are 
defined by analyzing the data separately for each two-day duty cycle. A table 
of coefficient values for each spherical harmonic representation is given along 
with global mean, gradients, degree variances, and contour plots. In addition, 
the entire data set is analyzed to define the monthly average radiation field. 


INTRODUCTION 

The importance of measuring and determining the Earth’s radiation budget 
(ERB) on global, zonal, and regional scales has been indicated by previous 
studies (refs. 1 and 2). Since the atmosphere is driven by radiation, a knowl- 
edge of the radiation budget is essential for understanding atmospheric processes 
(refs. 3 and 4). Much of our understanding of large-scale atmospheric processes 
was made possible by satellite measurements of Earth’s radiation budget and the 
ensuing analysis. 


Deconvolution by Parameter Estimation 

The ERB experiment aboard the Nimbus 6 spacecraft has provided valuable 
data from its wide field of view (WFOV) radiometers. Each data point is an 
integral of the irradiance from all points within the field of view of the 
WFOV sensor, a circular region on the Earth approximately 60° in diameter. 

F. B. House of Drexel University proposed that such data, being a convolution 
of the flux field at the top of the atmosphere, could be deconvoluted to enhance 
the resolution. Subsequently, the problem was solved by Smith and Green 
(refs. 5 to 7). Because of the statistical nature of the problem, a parameter 
estimation approach to the deconvolution problem was formulated. This tech- 
nique was simulated for both emitted and reflected radiation (wide and medium 
field of view). Parameter estimation verified the deconvolution concept and 
helped to define the propagation of errors (ref. 8). The eigenfunctions of the 
measurement operator were also shown (ref. 7) to be spherical harmonics for 
emitted radiation measurements, under certain conditions. Presented in this 
paper are discussions of the deconvolution estimation concept and results 
obtained by its application to the August 1975 Earth-emitted radiation data 
obtained with the WFOV ERB sensor aboard the Nimbus 6 satellite. 

The parameter estimation approach to deconvolution is followed in this 
paper. Basically, the radiation field at the top of the atmosphere is modeled 
in terms of a linear combination of base functions. The influence coefficients, 
or shape factors, are then computed. These influence coefficients describe the 



influence of model parameters on the radiation measurements at satellite alti- 
tude and enable the measurements to be expressed as a linear combination of 
the model parameters . This system of simultaneous measurement equations is 
solved by the Gauss-Markoff theorem. The parameter estimation approach is 
quite general and will accommodate various sensor models, spatially varying 
directional models for the radiation intensity, and elliptical orbital 
geometries . 

Under certain conditions, the eigenfunctions of the measurement operator 
are spherical harmonics. If these conditions are met, one can take advantage 
of this by defining the base functions, which model the radiation field, as 
spherical harmonics. The two conditions are: the sensor integrates incoming 

radiation over its field of view, with its directional response a function of 
nadir angle only; and the directional dependence of exiting radiation at any 
point at the top of the atmosphere is a function only of the zenith angle of 
the exiting ray. The first condition is satisfied by the ERB radiometer which 
is modeled by a perfectly black flat-plate sensor with a cosine nadir-angle 
response. The second condition states that the directional function is inde- 
pendent of latitude, longitude, and azimuth and is a function of the zenith 
angle only. That is, one directional function serves to describe the direc- 
tional dependence of exiting radiation for the entire globe. This is a reason- 
able assumption for longwave radiation. 

Not enough is currently known about the spatial variation of the direc- 
tional dependence of exiting radiation to model it effectively. Moreover, it 
has been shown in reference 8 that the gross features of the longwave radiation 
field are rather insensitive to this assumption. For these reasons, the base 
functions have been chosen as spherical harmonics. This approach has several 
advantages. The measurement operator integrates the incoming radiation over 
the sensor field of view. The influence coefficients are expressed in terms 
of this integration, which is normally performed by a numerical quadrature 
formula. This numerical integration requires considerable computational effort 
and introduces errors into the estimation process . Since spherical harmonics 
are eigenfunctions of the measurement operator, the problems associated with 
numerical integration are totally eliminated. The increase in accuracy, due 
to an exact quadrature formula, and the decrease in computational effort are 
felt to outweigh the consequences of modeling the directional response for 
longwave radiation as a function of the zenith angle only. Another advantage 
is that the radiation field is expressed in terms of spherical harmonics which 
are convenient for computing other descriptive quantities, such as the spatial 
spectrum of the radiation field. Thus, the parameter estimation technique is 
formulated in terms of spherical harmonics which simplify the deconvolution and 
facilitate its understanding. Since the Nimbus 6 orbit is circular (all mea- 
surements at the same altitude) and the measurement errors are modeled as inde- 
pendent random variables, the parameter estimation technique is further 
simplified. 


Background 

The Nimbus 6 ERB WFOV data for August 1975 have been analyzed by various 
investigators using different techniques. In the past, the data have been 
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analyzed on a monthly basis to yield an average radiation field for the month. 
This procedure, however, gives no information on how much the radiation field 
varies during the month. A monthly average is descriptive of the data only when 
the variations during the month are small. In this paper, the temporal varia- 
tion during the month is investigated by dividing the month into eight periods 
of four days each and analyzing the data for each period separately to produce 
a short time history of the Earth emitted radiation field. These results are 
then combined to obtain the monthly mean radiation field. An analysis of these 
data, reported by Smith et al. (ref. 9), defined the monthly mean radiation 
field by a global contour plot and the global mean. The data analysis tech- 
nique (height rectification) employed a geometric shape factor. - ! The same data 
were analyzed by Weaver (ref. 10) whose approach used the geometric shape fac- 
tor. It is worth noting the relationship between parameter estimation and 
height rectification whereby the radiation field at the top of the atmosphere 
is computed by multiplying the measurements by a single constant or shape fac- 
tor. The high order terms are attenuated by the WFOV measurement. Height 

rectification does not compensate for this attenuation. Parameter estimation, 
however, compensates for the high order attenuation and should lead to more 
accurate results. Another analysis of these data, by Green and Smith (ref. 11), 
defined the monthly mean radiation field by a 12th degree spherical harmonic 
contour plot . This analysis was based on the parameter estimation technique 
with the added assumption that the data can be averaged over 5° by 5° areas at 

the top of the atmosphere . This assumption reduces the computational effort 

but introduces error by smoothing the data. Green and Smith also present a 
15th degree spherical harmonic contour plot in reference 12. This analysis was 
actually extended to degree 36 to demonstrate the ill-posed nature of the decon- 
volution problem. Definition of the radiation field to degree 36 by the param- 
eter estimation approach, however, would require the inversion of a prohibi- 
tively large matrix. Therefore, the data were averaged over 5° by 5° areas and 
the radiation field was obtained by approximating the orthogonality of spherical 
harmonics. This approach is necessary to obtain estimates of the high frequency 
components of the radiation field, but it is not necessary or desirable for the 
low frequency components such as the global mean, pole-to-pole gradients, and 
pole-to-equator gradients. 


Objectives 

The purpose of this paper is to demonstrate the parameter estimation tech- 
nique with the Nimbus 6 ERB data. Since this approach yields the best (minimum 
variance) unbiased linear estimate of the radiation field within the modeling 
assumptions, it is felt to give more accurate results than previously reported. 
The analysis is described in detail and the results are presented in numerical 
form so that direct comparisons can be made. Previously reported results are 
inadequate for this purpose. 


1 Private communication . 
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SYMBOLS 


complex coefficients of spherical harmonics, W/m 2 
observation matrix 

vector of spherical harmonic coefficients, W/m 2 
real coefficients of spherical harmonics, W/m 2 
differential surface element 

directional function for emitted radiation, sr“1 

angular response of sensor 

satellite altitude above R, km 

radiant intensity, (W/m 2 )/sr 

number of measurements 

measurement operator (see eqs. (7) and (8)) 
measured flux at satellite altitude, W/m 2 
vector of measured fluxes 
degree of spherical harmonic expansion 

associated Legendre polynomial of degree n and order m 
radiative flux at top of atmosphere, W/m 2 

radius of surface approximating top of atmosphere (R = R e + 30), 

radius of Earth (6378.165 km) 

root mean square of measurement residual, km 

spherical harmonic representation of measured flux at satellite 
altitude, W/m 2 

complex spherical harmonic of degree n and order m 

real spherical harmonics of degree n and order m 

cone angle from satellite nadir to point on surface at top of 
atmosphere, deg 



aft cone angle from satellite nadir to horizon of surface at top of 

atmosphere , deg 

(3 clock angle from north about satellite nadir to point on surface 

at top of atmosphere, deg 

Y Earth central angle between satellite nadir and a point on surface 

at top of atmosphere, deg 

6 m Kronecker delta 

n 

e vector of measurement residuals, W/m 2 

e± ith measurement residual, W/m 2 

g azimuth angle of exiting ray, deg 

0 colatitude , deg 

X n nth eigenvalue of the measurement operator 

nth degree variance (see eq. (22)), W 2 /m^ 

4> longitude, deg 

\p zenith angle of exiting ray, deg 

0 ) solid angle, sr 

Superscripts : 

T transpose 

-1 inverse 

Subscript: 
s sensor 

A circumflex ( A ) over a symbol denotes an estimate. 


DATA 

The data used in this study to define the Earth’s longwave radiation field 
for August 1975 were obtained from the Earth radiation budget (ERB) instrument 
aboard the Nimbus 6 satellite. A description of the ERB instrument is given in 
reference 9. The data tapes were supplied by the National Oceanic and Atmo- 
spheric Administration. 
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Earth Radiation Budget Experiment 

The ERB instrument obtained both fixed wide-angle and scanning narrow-angle, 
measurements. The data considered here will be the fixed wide-angle irradiance 
at satellite altitude. One of the data channels recorded the total Earth irra- 
diance (0.2 to 50 pm) and another channel recorded the shortwave irradiance 
(0.2 to 3 • 8 ym) . The longwave contribution is the difference between these two 
measurements . The longwave irradiance data were then increased 1 1 percent as a 
calibration correction. A discussion of this correction is given in reference 9 
and accounts for a discrepancy between the ERB experiment fixed wide-angle and 
scanning narrow-angle measurements. The ERB instrument operated on a duty cycle 
of two days on and two days off. Irradiance measurements were taken every 
16 seconds. There were eight data periods during the month, but two of these did 
not have sufficient data to be analyzed separately. Approximately 50 000 mea- 
surements were gathered during August 1975, but not all of these data could be 
used in the analyses . 


Data Edit 

The data had to be edited to eliminate measurements which were unacceptable 
for various reasons. Some of the measurements were degraded by exposure to the 
Sun. In order for the wide-angle radiometer to view the Earth disk from horizon 
to horizon, the actual field of view was a few degrees larger than the Earth 
disk and included a thin annulus of deep space. When the satellite was near 
sunrise or sunset, the Sun came into view and contaminated the measurements. For 
this reason, measurements were omitted from the set of data when the Sun zenith 
at the subsatellite point was between 111.47° and 123.47°. Measurements were 
also omitted if they were greater than 280 W/m 2 or less than 50 W/m 2 . Similarly, 
a measurement which changed by more than 10 W/m 2 over a 1 6-sec time interval was 
considered unacceptable and was omitted from the data set. Another questionable 
set of data was obtained on August 16, 1975. Starting at 1000 hours GMT, the 
measurements gradually decreased in value during a full satellite revolution and 
then seemed to stabilize at values approximately 30 W/m 2 lower than expected. 

This occurred at the end of a data acquisition period. The measurements from 
the next data period did not exhibit this problem. The reason for this gradual 
decrease in value is not known. Thus, the measurements starting at 1000 hours 
GMT and continuing for three satellite revolutions, to the end of the data 
period, were omitted from the data set. The final edited data set contained 
45 942 measurements obtained during August 1975. 

During part of the orbit , the ERB radiometers viewed the Earth in dark- 
ness where the incident shortwave radiation is negligible. However, the direct 
irradiation of the sensor at spacecraft sunset caused a thermal transient in 
the measurements. As the field of view passed into darkness, the shortwave 
measurements rapidly decreased to about 9 W/m 2 in the first 2 min and continued 
to decrease during the remaining 30 min to a lower limit of 1 to 4 W/m 2 . No 
attempt was made to correct for this transient in the data. 
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ANALYSIS 


The edited ERB data were used to estimate the .radiation field at the top 
of the atmosphere. The radiation field was modeled as a linear combination of 
spherical harmonics and the parameters estimated with a least squares estimator. 


Radiation Model 

The top of the Earth-atmospheric system was approximated by a sphere of 
radius R. The intensity I of radiation leaving any point E on this surface 
is a function of colatitude 0, longitude <t>, and zenith angle ^ of the exit- 
ing ray. (See fig. 1.) The longwave radiative flux Q at point E is given by 


7T 

/"»£ =2it j|>=- 

Q(0,4>) = 1 \ 2 I C0 »4> ) cos sin dtp <3C (1) 

J C=0 J\Ij=0 


where £ is the azimuth of the exiting ray. The intensity was modeled as 


I(0,<i>,\|O = 1 Q(0 ,4>) F(4>) 

7T 


( 2 ) 


where F(\p) is a directional function and is dependent upon only the zenith 
of the exiting ray. Notice that this function defines the directional charac- 
teristics of the exiting radiation at every point on the surface. The strength 
of the radiation at a given point, however, varies with colatitude and longitude 
and is defined by Q(0,<f)). Substituting equation (2) into (1) yields the 
normalization condition for F(^), that is, 
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F (\p ) cos y\) sin d\p = 1 


(3) 


If the surface radiation is Lambertian, that is, if radiation leaves the surface 
with equal intensity in all directions, then F(\p) = 1. If this is not true, 
then the directional characteristics can be defined by a limb darkening function, 
an example of which is given by 
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1.074 exptO. 106(1 - sec ip)] 

FOIO 1.074 exp [-0 .056 + 0.05(1 - sec \p)] 


(0° ^ < 60°) 

( 60 ° ^ < 90 °) ( 4 ) 

( 4 > > 90 °) 


This function results from averaging the family of limb darkening functions 
in reference 13. 


Measurement Model 

Now consider the radiation sensor at satellite altitude h, as shown in 
figure 1 . The radiation from the Earth incident on the sensor at a colati- 
tude 9 S and a longitude 4> s is 


m(6 s ,<t> s ) 


I 

FOV 


1(9, <J>»40 g(a) du 


(5) 


where w is the solid angle at the satellite subtended by the surface element, 
a is the cone angle at the satellite from the local vertical to the surface 
element, and the integration is over the field of view (FOV) of the sensor. 

The function g(a) is the angular response of the sensor to incoming radiation 
and is modeled as a perfectly black flat-plate sensor, or 


g(a) = cos a 


( 6 ) 


Substituting (6) and (2) into (5) yields 


ro(9 s ,<f> s ) = - \ Q(9,tj>) F(\|>) cos a dio 

IT O 

FOV 


or 


1 pB=2ir r»a=a h 

m(9 s ,4>s) = - \ \ Q(9,(j>) F (\p) cos a sin a da dg (7) 

v J g=0 J a=0 
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where is the cone angle from the satellite nadir to the horizon. It 

is convenient to express this relationship between the measured radia- 
tion m(0 s ,<t> s ) and the unknown flux Q(0,4>) as 


m(0 s ,<t> s ) = L[Q(0,4>)] 


( 8 ) 


where L denotes the linear integral measurement operator of equation (7). 
Smith and Green (ref. 7) have shown that the eigenfunctions of this linear 
operator are spherical harmonics, that is, 

L [z” (e ,4>)] = X n Z®(0 S ,4> S ) (9) 


where Z™(0 S ,4> S ) is a spherical harmonic of order m and degree n evaluated 
at the subsatellite point. The associated eigenvalue X n is given by 


C a=a h o 

X n = 2 \ P (cos y) F(4>) cos a sin a da (10) 

^a=0 n 

where P^(cos y) denotes the Legendre polynomial of degree n as a function 
of the Earth central angle y . (See fig. 1.) 


Deconvolution 

Because spherical harmonics are eigenfunctions of the measurement oper- 
ator, it is convenient to let the flux at the surface be represented by a 
truncated series of spherical harmonics, that is, 


N n 

0(0,4.) = ^ ^ b n z n (0 ’ <1)) (11) 

n=0 m=-n 


Also, let the measurements be represented as 


m(0 ,4>) 


I i 

n=0 m=-n 


( 12 ) 
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Substituting equations (11) and (12) into (8) gives 


N n 


Yj X Z n (0s,<t)s) = L 


n=0 m=-n 


but from equation (9) 


N n 




n=0 m=-n 


N n 


N n 


£ £ a” Z”(0 s »4)s) = X Z Xn b n Z n (9s ^ s) 


n=0 m=-n 


n=0 m=-n 


Equating coefficients of like terms gives 


or 


m , .m 
a = X b 
n n n 


m 


m 


b = — 
n \ 
A n 


(13) 


Thus, the flux at the surface is 


N n 


0(6,4,) = ^ ^ b n Z n (0 ’ ( t ,) 


n=0 ro=-n 


or 


N n a m 


Q(0 


■ 2 I r n Z > A) 


(14) 


n=0 m=-n 
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It is seen that the flux depends on the coefficients of the measurement repre- 


sentation a n and the eigenvalues of the measurement operator X n . 


Spherical Harmonics 

For the purposes of this analysis, a real formulation of the spherical 
harmonics was used. Thus, if X(0,<t>) is a real function over a sphere, then 
a representation of X(0, <f>) is 


N n 

X(e,<i>) = 2 j [ C n + S n Z s n (0 ’ <|)) ] (15) 

n=0 m=0 


where Z™ C0 »<J> ) are spherical harmonics given by 


O 0 ^ = 


(2n + 1 ) (n - m) 




(n + m) ! 


1/2 


cos mch P m (cos 0) (16a) 

n 


Z m (0,cj>) = 
sn 


(2n + 1 ) (n - m) 


! (* - 0 


(n + m) ! 


1/2 


sin m<{) P™(cos 0) ( 1 6b ) 


and P^Ccos 0) are associated Legendre polynomials. These polynomials are 
defined by the following recursive equation: 

P 0 (cos 0) = 1 
0 

0 

P.j(eos 0) = cos 0 
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P . ( cos 0 ) 
n+1 



0 0 

+ 1) cos 0 P (cos 0) - nP .(cos 0) 
n n-1 


>] 


(n = 1, 2, 


P m (cos 0) = (2m - 1) sin 0 P m- !](cos 0) 
m m-1 


= (2m + 1) cos 0 P m (cos 0) 
m 

(2n +1) cos 0 P m (cos 0) - (n + m) P m .(cos 0) 
n n-1 

n - m + 1 

( m = 1 , 2, . . 

n = m + 1, m + 2, . . 

It follows that these spherical harmonics are orthogonal according to 

f z“ (0,4) zt . (0 ,4>) do = f z“ (0,6) Z* (0,4) da = im 

xj cn cj xj sn sj 

Sphere Sphere 


P™. .(cos 0) 
m+1 

P m . ( cos 0 ) 


except for 


5 z sn (e '' t,) Z S 0 n <e - 4) da = 0 
Sphere 

Also , 

J z“ n (0,4) Z sj (9 » <1)) da = 0 

Sphere 



fi 1 
n m 

(17a) 


(17b) 


(17c) 
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With these definitions, the flux at the surface (eq. ( 1 -4 ) ) becomes 


Q(e,4>) 


N n 

1 1 =R 


n=0 m=0 



where C m and S m 
n n 


are the coefficients of the measurement representation. 


Estimation of Coefficients 

Given a discrete set of measurements , it now remains to determine the 
coefficients of the measurement representation. These coefficients are esti- 
mated with a least squares fit to the measurements . Each discrete measure- 
ment m(0j_,({)i) is expressed as a linear combination of the coefficients; thus, 
the ith measurement is given by 


N n 

®(ei,4>i) = ^ £ [ C n C (0i ^ i) + S n + ei (18) 

n=0 m=0 


where is the measurement residual, or the difference between the actual 

measurement and its representation. This equation and the estimation equations 
to follow can be greatly simplified by ordering the terms such that the coef- 
ficient vector is defined as 


C = 






(19) 


which consists of (N + I) 2 terms. Note that S° is not contained in the 

n 

ordering. From equation (16b), it can be seen that the spherical harmonic 

m 0 

z“ (0,<t>) = 0 when m = 0. Hence, its coefficient S is meaningless. Also 
sn n 

define a measurement vector as 


M = [10(0! ,<|>-|), m(02,<l>2)* 


m (0k» T 
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a residual vector as 


e = (e-|, e 2 , 


e k ) 


and an observation matrix as 


Z q 0 (e 1 *<(> 1 ) z° 1 (0 1 »4>i) 


Z cO (02,<i)2) Z cl (02 ’ <t ’ 2) 


B = 


Z sN <e '" l> ' ) 

Z 3N <e 2-*2> 


• • • ^sfj^kj^k) 

_ 


where k is the number of measurements . The entire set of measurements given 
by equation (18) is now conveniently expressed in matrix notation by 


M = BC + e 


( 20 ) 


This is an overdetermined system of simultaneous equations. The solution that 
minimizes the sum of squares of the residuals e T e is the least squares solu- 
tion and is given by (ref. 14) 


6 = (B t b)-‘> B t M 


( 21 ) 


RESULTS 

The ERB data were analyzed to demonstrate and evaluate the concept of 
parameter estimation and to define the variation of the longwave radiation 
field during the month of August 1975. The spatial variation is given by 
the spherical harmonic functions in colatitude and longitude. The temporal 
variation was obtained by analyzing the data separately for two day periods. 
This is a natural division of time since the ERB instrument operated on a 
duty cycle of two days on and two days off. During the month of August 1975, 
there were eight duty cycles, two of which had an insufficient number of 
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measurements to be analyzed separately. Thus, the data were analyzed sepa- 
rately for six duty cycles of two days each, resulting in a short time history 
of Earth's longwave radiation. In addition, the entire data set for all periods 
was combined and analyzed to determine a monthly average radiation field. 


Degree Variance 

The data obtained during the first duty cycle (August 2-4) were analyzed, 
and the radiation field at the surface was defined by a 12th degree and 
12th. order spherical harmonic expansion where the coefficients were least square 
estimates (eq. (21)). These coefficients define the spatial spectrum of the 
radiation field in terms of the degree variance. 

The degree variance is defined as the average value of the solution 
squared for a given degree. If X(0,4>) is a spherical harmonic representation 
(eq. (15)), the degree variance is given by 


s 


2 Sphere n 


X“(0,<t>) da 


n 


X 


da 


Sphere 


where 


( 22 ) 


n 

x n ( M> = £ [c" z“ n <e,4>) . s” z” n <e,*)] (23) 

m=0 

Substitution of equation (23) into equation (22) along with the orthogonality 
conditions (eq. (17)) yields 


°n ■ l feT * «f 
m=0 L 


where is defined as zero. The square root of degree variance is called 

degree dispersion (denoted by o n ) and is often used to denote the variation. 
The degree variance for all n is analogous to the amplitude spectrum in ordi- 
nary harmonic analysis and gives an indication of the power in each frequency <, 
The degree variance for the first duty cycle is shown in figure 2 at both sat- 
ellite altitude and at the surface. The relationship between the two degree 
variances follows from equation (13) and is given by 


15 



2 


a 


n 


(surface) 


2 

ct (altitude) 
n 


X 


2 

n 


where \ n are the eigenvalues defined by equation (10). Eigenvalues are pre- 
sented in table I for a Lambertian directional model and for the limb darkening 
model defined by equation (4). The sequence of eigenvalues is monotonically 
decreasing and results from the smoothing effect of the measurement operator. 
Since the amplitude of the high frequency components of the radiation field are 
reduced by the measurement operator, the inverse process must amplify these com 
ponents. This can be seen by the divergence in the two degree variances as the 
degree number increases. Representing the radiation field by increasingly 
higher degree spherical harmonics will eventually lead to a meaningless repre- 
sentation at the surface due to this divergence. Variations of the radiation 
field and errors in the measurement produce high frequency components at satel- 
lite altitude which are greatly amplified when converted to the surface by the 
factor 1/Xn* Thus, there is some degree beyond which variability effects and 
measurement errors dominate. The representation for the first duty cycle, how- 
ever, seems to be well behaved at least to degree five. Therefore, the radia- 
tion field for this study was modeled by a fifth degree and fifth order expan- 
sion in spherical harmonics. 


Estimated Coefficients 

The coefficients of the radiation field at satellite altitude are pre- 
sented in table II for the first six duty cycles in August. Each radiation 
field is represented by the 36 coefficients of the spherical harmonic expansion 
The ordering corresponds to the scheme defined in equation (19). The last 
column is the monthly average radiation field and covers the period August 2 
through August 31* The data set for this radiation field contained JJ5 9^2 mea- 
surements and incorporated all of the measurements from each of the six duty 
cycles plus 1978 measurements from the eighth duty cycle. This duty cycle con- 
tained good data for about three orbital revolutions which are beneficial to 
the monthly average but not sufficient to be analyzed separately. The seventh 
duty cycle contained no data. Also presented in table II is the RMS residual 
between the measurement and the spherical harmonic representation, that is, 


RMS 



K 

I [ m(e i j <t>i) ~ x(ei,<H)] 

i=1 
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where k is the number of measurements. According to equations (18) and (20), 
the RMS is expressed in matrix notation as 


RMS 


= [-I eT£ ] 1/2= [>- B£)T 


-l1/2 


(M - BC) 


Expanding this expression and incorporating equation (21), yields 


RMS 2 r -[m t M - C T B t m] 


The RMS is a measure of the amount of total variation M T M that has not been 
accounted for by the spherical harmonic representation. It is, therefore, an 
indication of how well the measurements are represented by X(0,<t>). 

The spherical harmonic representation X(0,4>) of the measurements at satel- 
lite altitude is deconvoluted to the surface by dividing by the eigenvalues of 
the measurement operator (eq. (13)). This transformation is based on the 
assumed directional function of the radiation at the surface. The radiation 
field at the surface for Lambertian radiation is presented in table III, and 
the degree dispersions are presented in table IV. Some of these coefficients 
have obvious physical meanings. For example, the global mean flux is given by 


Qglobe 


5 

Sphere 


Q(9,4>) do 



Sphere 



Sphere 


IIES Z cn (e ’+ ) * S n Z sn (9 '*>] 
n=0 m=0 


da 
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All spherical harmonics integrate to zero over the surface except Z c q(6,<|>) o 1 
which integrates to 4n . Thus, the first coefficient is the global mean 

flux at the surface. In addition, the coefficient of the first zonal har- 


monic C;j is a measure of the pole-to-pole gradient, and the coefficient of 
the second zonal harmonic is a measure of the equator-to-pole gradient . 


Interpretation of Coefficients 


The radiation field for each duty cycle is depicted by a contour plot in 
figures 3 to 8, and the monthly average radiation field is depicted in figure 9- 
These representations should not be interpreted as final since they are only 
fifth degree representations of the radiation field. The data should support a 
higher degree representation (fig. 2). Nevertheless, a fifth degree represen- 
tation is sufficient to demonstrate the parameter estimation technique and 
define the variation of the major features during one month. From the contour 
plots, it can be seen that the positions of the major highs and lows are rather 
stable during the month. For this reason, the monthly average plot (fig. 9) is 
representative of the major features of the radiation field throughout the month 
It can also be seen from the contour plots and from the table of coefficients 
that the radiation fields are zonal in character. This can also be seen by com- 
puting the present contribution of the zonal coefficient to the degree variance, 


that is, 



100 $. 


For the monthly average radiation field, the percent- 


ages for increasing degree numbers are as follows: 100, 85, 90, 76, 70, and 17. 

For degree zero , the entire degree variance is composed of the zonal contribu- 
tion. The zonal component is also the dominant influence for degrees one to 
four. Degree five shows more of a balance between the zonal and nonzonal terms. 
Since the degree variance generally decreases with increasing degree (at least 
through degree five for this case), the conclusion is reached that the radiation 
fields are predominately zonal in character. 


Insight into the nature of the radiation fields can be obtained by order- 
ing the 36 coefficients according to their importance. This is accomplished by 
computing the increase in the RMS value when a specific term is omitted, that is 
the importance of a term is equated to the increase in - the RMS when a 35 coef- 
ficient representation is determined. This procedure revealed that the five 
most important coefficients for the first duty cycle were the first five zonal 
harmonic coefficients. The 6 zonal harmonic coefficient numbered 10th in impor- 
tance. Generally, a size ordering of the coefficients follows closely the 
importance ordering of the coefficients. If the position of the measurements 
were such that the B matrix in equation (21) was diagonal, then the two 
orderings would be identical. This is not true for the Nimbus 6 data since the 
measurement positions correspond to satellite positions. Nevertheless, for the 
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data under consideration, the two orderings differed in the first 10 positions 
only by a rearrangement of the 7th, 8th, and 9th terms. 


An alternative to using a complete set of spherical harmonics (i.e., using 
all the terms in the sequence of eq. (19) for a given N) is to use a selected 
set of harmonics. A subset of the most important harmonics would be appropriate. 
Consider the first duty cycle. Elimination of all terms except the global mean 

flux Cq corresponds to a one term representation with an RMS value of 29.916. 


The total variation is defined as the square of this quantity. From table II 
the RMS value for a 36 term representation is 9.564. Thus, a 36 term solution 


accounts for 


(29.916)2 - (9.564)2 


100/6 s 90% 


of the total variation. 


Like- 


L (29.916)2 J 

wise, a solution based on the 20 most important terms accounts for 88 percent 
of the total variation and a 10 term solution accounts for 84 percent. Hence, 
the 10 most important terms account for 84 percent of the total variation, the 
next 26 terms in importance account for 6 percent, and the 16 least important 
terms account for only 2 percent. 


Importance of Directional Model 


When the radiation field at satellite altitude is deconvoluted to the 
surface by equation (13), a directional model must be defined through the 
eigenvalues. Sets of eigenvalues are presented in table I for a Lambertian 
model and a limb darkening directional model. The surface radiation fields 
in table III are based on the Lambertian model. The difference in the sur- 
face radiation field resulting from the two different directional models is 
very small especially for the low frequency terms. The global mean is defined 


by the constant term given by 


0 1 

C n ( surface) = — 
u Xn 


^(altitude) . 


But, Xq is the 


same for both directional models and, in fact, is the same for any directional 
model defined by the normalizing condition (eq. (3)). Thus, the global mean 
at the surface is not a function of the directional model. For a fifth degree 
representation, only the first six eigenvalues are needed. These eigenvalues 
for the two models differ only slightly. Moreover, most of the power is in 
the lower frequencies. Hence, the fifth degree representation is not a strong 
function of the directional models. 


Time History of Coefficients 

A short history of the most important coefficients is presented in fig- 
ure 10. The change in the global mean estimate from its least value to 

its largest value is only 2.0 W/m2. Its greatest departure from the monthly 
mean is 1.3 W/m2. It is seen that the time variations in the other coefficients 
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are also small. Because they are so small, it is not clear whether they are 
due to real fluctuations of the atmosphere, sampling variations, or some other 
cause . 


CONCLUDING REMARKS 

The Nimbus 6 Earth radiation budget data for the month of August 1975 has 
been analyzed by the parameter estimation technique to define the variation of 
the longwave radiation field at the top of the atmosphere. Since the eigen- 
functions of the measurement operator are spherical harmonics, they provided a 
convenient representation of the radiation field. Moreover, the coefficient of 
the first zonal harmonic is the global mean flux, while the coefficients of the 
second and third zonal harmonics are measures of the pole-to-pole gradient and 
the equator-to-pole gradient, respectively. It was shown that the zonal har- 
monics are the most important terms in the representation of the radiation 
field; that is, the radiation field at the top of the atmosphere is character- 
ized by axial symmetry. 

The radiation field was represented as a fifth degree and fifth order 
expansion in spherical harmonics. Atmospheric variations and random measure- 
ment errors cause the representation to diverge in the high frequency compo- 
nents . The degree variance, however, showed that the solution was well behaved 
through degree five which corresponds to representing the solution in terms of 
36 coefficients. Even though a higher degree representation can be supported 
by the data, a fifth degree solution is sufficient to demonstrate the parameter 
estimation technique and define the variation of the major features during the 
month. If the coefficients are ordered according to their importance in reduc- 
ing the root mean square of the measurement residual, the radiation field can 
be expressed in terms of a subset of those coefficients. When all 36 coeffi- 
cients are used, they account for 90 percent of the total variation. A repre- 
sentation based on the 20 most important coefficients accounts for 88 percent 
of the total variation, and a 10 coefficient representation accounts for 84 per- 
cent. Thus, the 10 most important coefficients account for 84 percent of the 
total variation, the next 26 coefficients in importance account for 6 percent, 
and the 16 least important coefficients account for only 2 percent. 

The data were analyzed for two-day periods of time resulting in a short 
time history. The coefficients corresponding to the low frequency zonal har- 
monics are fairly constant with time. For example, the change in the global 
mean estimate from its least value to its largest value was only 2.0 W/m^, 
or less than 1 percent . Its greatest departure from the monthly mean was 
1.3 W/m^. Because these changes are small, it is not clear whether they are 
due to real fluctuations of the atmosphere, sampling variations, or some other 
cause. 

When the radiation field at satellite altitude is deconvoluted to the top 
of the atmosphere, the directional characteristics of the radiation field must 
be defined and are expressed through the eigenvalues of the measurement oper- 
ator. Examination of these eigenvalues shows that the representation to fifth 
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degree is not strongly affected by these characteristics. In fact, the global 
mean is seen to be independent of the directional model. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
September 19, 1978 
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TABLE I.- EIGENVALUES OF MEASUREMENT OPERATOR 

|"R = R e + 30.000 km = 6408.165 km;l 
[_ h = 1 100 km - 30 km = 1070 km J 


n 


for - 

Lambertian 

model 

Limb darkening 
model 

0 

0.7343 

0.7343 

1 

.7217 

.7232 

2 

.6975 

.7014 

3 

.6632 

.6704 

4 

.6208 

.6317 

5 

.5726 

. 5873 

6 

.5214 

.5393 

7 

.4693 

.4899 

8 

.4185 

.4408 

9 

• 3707 

• 3936 

10 

• 3267 

.3494 

11 

.2874 

.3091 

12 

.2526 

.2728 
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TABLE II.- COEFFICIENTS OF SPHERICAL HARMONIC REPRESENTATION AT SATELLITE ALTITUDE 







Coefficients for duty cycle - 



n 

m 

Coefficient 

1 

2 

3 

4 

5 

6 

Monthly 

average 

0 

0 

C 

173.024 

172.430 

172.319 

173.212 

173.826 

173.106 

172.860 

1 

0 

C 

8.911 

8.058 

9.060 

8.550 

8.586 

7.696 

8.535 

1 

1 

C 

3.624 

3.662 

2.598 

2.663 

2.828 

1 .766 

2.689 



S 

-2.739 

-1.768 

-2.759 

-3.197 

-1.653 

-1.307 

-2.358 

2 

0 

C 

-15.906 

-16.088 

-14.980 

-16.782 

-16.651 

-16.627 

-16.074 

2 

1 

C 

2.229 

4.066 

3.904 

2.042 

2.581 

2.602 

2.777 



S 

-.731 

-1.670 

-1.763 

-.800 

-.832 

-1 .518 

-1 .203 

2 

2 

C 

2.815 

4.044 

3.696 

3.693 

2.896. 

3.589 

3.477 



S 

3.606 

2.715 

3.842 

1.913 

2.162 

1 .223 

2.619 

3 

0 

c 

8.413 

9.196 

9.227 

8.009 

8.338 

7.603 

8.473 

3 

1 

c 

.068 

.469 

1.133 

1.233 

1.533 

.395 

.853 



s 

3.461 

3.353 

3.663 

3.502 

3.571 

3.132 

3.392 

3 

2 

c 

2.454 

2.343 

2.043 

2.125 

1.401 

1.633 

1 .997 



s 

2.223 

3.545 

1 .894 

2.130 

3.256 

1.874 

2.358 

3 

3 

c 

1 .620 

.514 

.975 

1.345 

-.052 

-.678 

.704 



s 

-1 .586 

.861 

-.879 

-.574 

-.350 

-1.357 

-.610 

4 

0 

c 

-6.177 

-5.883 

-6.054 

-4.850 

-6.422 

-5.971 

-5.984 

4 

1 

c 

-.918 

-1.143 

-.537 

-.270 

-1 .082 

-.248 

-.599 



s 

3.825 

2.875 

3.212 

2.523 

1.797 

1.178 

2.637 

4 

2 

c 

-1 .534 

-.213 

-2.121 

-2.538 

-2.245 

-.540 

-1.537 



s 

.143 

1.093 

1.143 

-.024 

-.073 

.302 

.394 

4 

3 

c 

-1.379 

-2.007 

-1.923 

-.302 

.005 

-.984 

-1 .008 



s 

-.145 

.370 

-.152 

-.506 

-2.420 

.348 

-.278 

4 

4 

c 

-1 .069 

-1 .426 

-2.020 

-1 .912 

-2.050 

-2.892 

-1.992 



s 

.402 

-.726 

.901 

.589 

-1 .078 

.043 

.115 

5 

0 

c 

-2.866 

-1.936 

-2.530 

-1 .667 

-1 .552 

-1.104 

-2.053 

5 

1 

c 

-1.145 

-1 .456 

-2.258 

-2.379 

-2.772 

-1 .708 

-1 .883 



s 

-2.081 

-3.241 

-2.875 

-3.209 

-3.015 

-1 .692 

-2.645 

5 

2 

c 

-1 .325 

-.725 

-1 .221 

-2.209 

-1 .963 

-2.114 

-1 .555 



s 

-.332 

-.557 

-1.583 

-.630 

-.602 

.607 

-.506 

5 

3 

c 

-.074 

-.272 

.901 

-1.128 

.554 

-.253 

-.079 



s 

2.429 

1.625 

1 .289 

.638 

1.785 

2.175 

1.695 

5 

4 

c 

-1.131 

-1.077 

-1 .726 

.249 

-.233 

-2.198 

-1.090 



s 

-.075 

-1 .681 

-.341 

.153 

1 .267 

-.050 

-.209 

5 

5 ' 

c. 

-1 .355 

-2.068 

-1.387 

-2.310 

-1 .540 

.952 

-1.336 



s 

J 

-1.040 

.274 

-1 .412 

-1.394 

-1.108 

-.396 

-.886 










First 

day 



August 2 

August 6 

August 10 

August 14 

August 19 

August 22 

August 2 

Last 

day . . 


August 4 

August 8 

August 12 

August 16 

August 20 

August 24 

August 31 

Number of measurements . . 

8796 

6465 

8696 

7625 

5469 

6913 

45 942 

RMS residual 


9.564 

9.108 

9.467 

9.135 

8.987 

9.091 

10.016 
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TABLE III.- COEFFICIENTS OF SPHERICAL HARMONIC REPRESENTATION 

AT TOP OF ATMOSPHERE 


n m Coefficient 


Coefficients for duty cycle - 


235.632 

C\J 

C\J 

00 

^T 

on 

CO 

234.671 

235.887 

236.723 

235.744 

12.347 

11.166 

12.554 

11.847 

11.897 

10.663 

5.021 

5.074 

3.600 

3.690 

3.918 

2.447 

- 3.795 

- 2.450 

- 3.824 

- 4.430 

- 2.290 

- 1.810 

- 22.805 

- 23.065 

- 21 .476 

- 24.060 

- 23.872 

- 23.838 

3.196 

5.830 

5.597 

2.927 

3.700 

3-730 

- 1.048 

- 2.394 

- 2.528 

- 1.146 

- 1.193 

- 2.177 

4.036 

5.798 

5.299 

5.295 

4.153 

5.145 

5.171 

3.892 

5.509 

2.743 

3.100 

1.753 

12.685 

13-866 

13.913 

12.077 

12.572 

11.464 

.102 

.707 

1 .708 

1 .859 

2-311 

.595 

5.219 

5.056 

5.524 

5.281 

5.385 

4.723 

3.700 

3.533 

3.080 

3.205 

2.113 

2.462 

3.352 

5.345 

2.856 

3.212 

4.910 

2.826 

2.442 

.775 

1.470 

2.028 

-.078 

- 1 .022 

- 2.391 

1.299 

- 1.326 

-.865 

-.527 

- 2.046 

- 9.951 

- 9.476 

- 9.751 

- 7.813 

- 10.344 

- 9.618 

- 1.479 

- 1.841 

-.865 

-.435 

- 1.743 

-.399 

6.161 

4.632 

5.174 

4.065 

2.894 

1.898 

- 2.471 

-.342 

- 3-416 

- 4.089 

- 3-616 

-.870 

.231 

1.760 

1 .840 

-.039 

-.118 

.487 

- 2.221 

- 3.233 

- 3.098 

-.487 

.009 

- 1 .585 

-.234 

.596 

-.244 

-.815 

- 3.298 

.560 

- 1.722 

- 2.297 

- 3.253 

- 3.080 

- 3.303 

- 4.659 

.648 

- 1.169 

1 .452 

.949 

- 1.737 

.069 

- 5.005 

- 3.381 

- 4.418 

- 2.911 

- 2.710 

- 1.929 

- 2.000 

- 2.542 

- 3.944 

- 4.155 

- 4.840 

- 2.984 

- 3.634 

- 5.661 

- 5.020 

- 5.604 

- 5.265 

- 2.956 

- 2.314 

- 1.266 

- 2.133 

- 3.858 

- 3-428 

- 3.692 

-.580 

-.973 

- 2.765 

- 1.099 

- 1.051 

1 .060 

-.129 

-.475 

1.574 

- 1 .970 

.968 

-.441 

4.242 

2.837 

2.250 

1.114 

3-117 

3.799 

- 1.976 

- 1.881 

- 3.015 

.434 

-.407 

- 3.839 

-.131 

- 2.936 

-.595 

.267 

2.212 

-.088 

- 2.366 

- 3.612 

- 2.423 

- 4.034 

- 2.690 

1.662 

- 1 .816 

.478 

- 2.465 

- 2.435 

- 1.936 

-.692 


- 9.639 

-.965 

4 . 21)8 

- 2.477 

.635 

- 1.623 

-.442 

- 3.209 

.185 

- 3.585 
- 3.288 
- 4.619 
- 2.716 
-.884 
-.139 
2.960 
- 1 .904 
-.365 
- 2.334 
- 1.547 








Figure 1.- Earth-satellite geometry. 
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Figure 3 •- Map of longwave radiation at top of atmosphere for August 2-4, 1975. 
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Figure *J.- Map of longwave radiation at top of atmosphere for August 6-8, 1975. 
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Figure 6.- Map of longwave radiation at top of atmosphere for August 1*1-16, 1975. 
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Figure 7.- Map of longwave radiation at top of atmosphere for August 19-20, 1975. 
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Figure 8.- Map of longwave radiation at top of atmosphere for August 22-24, 1975. 
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Figure 10.- Time variation of zonal coefficients for longwave radiation at top of 

atmosphere for August 1975. 
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